The quantum scattering time and its implications on scattering sources in graphene 
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We determine the quantum scattering time in six graphene samples with mobility of 4,400 
< /X < 17,000 cmVVs over a wide range of carrier density (1.2 < n < 6xl0^^/cm^). derived from 
Shubnikov-de Haas oscillations ranges ~25-74 fs, corresponding to a single-particle level broadening 
of 4.5-13 meV. The ratio of the transport to quantum scattering time rt/rq spans 1.5-5.1 in these 
samples, which can be quantitatively understood combining scattering from short-ranged centers 
and charged impurities located within 2 nm of the graphene sheet. Our results suggest that charges 
residing on the Si02 surface play a dominant role in limiting carrier mobility in current samples. 

PACS numbers: 73.63.-b, 73.21.-b, 73.43.-f, 72.15. Lh 



Understanding and eliminating extrinsic scattering 
sources in graphene is critical to the advancement of its 
fundamental study and technological applications. De- 
spite many theoretical and experimental investigations 
into possible candidates, including charged impurities 
(CI), adsorbates, substrate corrugations, and ripples, 
contradictory observations remain and a clea r picture has 
yet to cmcrgc .'^^''^'^'^''^'^'^'^"'"'^^'^^'^'^'^^'^ 

To date, most experimental studies have focused 
on probing the carrier mobility /i, or equivalentl y the 
transport scattering time Tt = m*/x/e.'^'^° ^"^ I 3|i4 | i5 | 
Another important parameter in two-dimensional (2D) 
transport, the quantum scattering time r^, has not been 
well studied. Tq characterizes the momentum relaxation 
of a quasi-particle and relates to its quantum level broad- 
ening r through r — hj^Tq. 

Quantitatively, the difference between Tq and in 
graphene is shown in the following equationsl^^ 

1 = / Q(6')(l-Hcos6')d6i (1) 

1 r 

Q(6i)(l + cos6i)(l-cos6i)d6i. 

Tt Jo 

Here, 9 is the scattering angl e and Q{0) depends on 
specific scattering mechanisms ! ^ ^ 1 While small-angle 
events weigh heavily towards Tq, Tt is mostly affected by 
right angle scatterings. Measurement of Tf/Tq has proven 
to be a powerful diagnostic tool in revealing complex 
scattering scen ario s in conventional 2D electron gases 
(2DEGs).'^^ 2-fll2U por example, short-ranged scattering 
sources give rise to rt/r^ ~ 1 while charged impurities 
far away from a 2DEG lead to predominately small-angle 
events, resulting in large Tt/Tq. The former has been ob- 
served in silicon inversion layers and the latter character- 
izes modulation doped GaAs 2DEGs.^® ^° Despite its 
demonstrated importance, the study of Tq in graphene 
has been scant. Existing data are largely obtained from 
the linewidth of cyclotron resonance at low densities.'22l 
A systematic comparison between Tt and Tq has not been 
made. 

In this work, we report a comprehensive study of Tq 
in six graphene samples over a wide range of carrier 



densities 1.2 < n < 6 x 10^^/cm^ and mobility 4,400 
< /i < 17,000 cm^/Vs. Tq is obtained from Shubnikov-de 
Haas (SdH) oscillations and ranges approximately 25-74 
fs in these samples, corresponding to F = 4.5-13 meV. 
The n-dependence of tj, Tq, and their ratio Tt/Tq can all 
be ex plain ed by a self-consistent Boltzmann transport 
theorjEE^ using three parameters: the charged impurity 
density riimp, the impurity-graphene distance z, and the 
resistivity from short-ranged scatterers p short- Our re- 
sults indicate that the mobility in current graphene-on- 
Si02 samples is limited by scattering from charges re- 
siding within 2 nm of the graphene sheet. We speculate 
that charges present at the graphene/Si02 interface are 
the major sources of scattering. 

Single-layer graphene sheets are mechanically exfoli- 
ated onto 290 nm Si02/doped Si substrates and identi- 
fied optically. Rectangular pieces are processed into Hall 
bar devices using standard e-beam lithography followed 
by metal deposition (Fig. [l] inset). The fabrication de- 
tails are given in Ref.^'^. Representative data from four 
samples (denoted as samples A-D) are presented in de- 
tails. 

Transport experiments are performed in a pumped He^ 
cryostat with a base temperature of 1.4 K and equipped 
with a 9 T magnet. Standard lock-in techniques are used 
with an excitation current of 50-200 nA. The doped Si 
substrates serve as back gate electrodes, to which a volt- 
age (Vg) is applied to tune the carrier density and hence 
the conductance of graphene. We extract carrier density 
from SdH oscillations and obtain a gating efficiency of a 
= dn/dVg = 7xlOi°/cm2V of the backgate. All mea- 
surements are taken at T < 10 K to eliminate electron- 
phonon scattering. 

Figure [Tj shows the conductivity a vs Vg taken on sam- 
ples A-D. o'(Vg) varies linearly with Vg for the whole 
density range in sample A and exhibits sublinear cr(Vg) 
at high Vg for samples B-D, where the "bowing" is most 
pronounced in sample C. The qualitative fe atures of t hese 
traces resemble results reported previously,'^^!!!!!!!!!! ^^f^ 
have been explained by a self-consistent Boltzmann the- 
ory combining scattering from long and short-ranged 
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FIG. 1: (Color online) o-(Vg) of samples A (black), B (red), 
C (green), and D (blue) below 10 K. Inset: Optical image of 
sample A. Edge of the piece outlined. Scale bar is 10 ^im. 
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TABLE I: fiFE, Pshort, Tq, Tt/Tq, z, and Tiimp for samples A-D. 
Tq and Tt/rq are given for n ~ 3 x 10^^/cm^. The uncertainty 
in z is 1-2 A for all samples. Data in parenthesis are after the 
corrections of density inhomogeneity. 



sources 
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a 



mm 



Pshort', f^lo 



ne/ipE + CTr 



(2) 



In this framework, Pshort denotes a constant contribu- 
tion to resistivity from short-ranged scattering sources 
such as defects or neutral adsorbates. CI are thought to 
give rise to a linear (T(Vg), which implies a constant field 
effect mobility /xpE and consequently Tf cx ^/n. Equa- 
tions [2] produce excellent fittings to the 'T(Vg) data of all 
our samples. The resulting /ipE and Pshort span 4,400- 
17,000 cmVVs and 40-165 n respectively (Table |l|, cov- 
ering much of the variations reported in the lit eratu re. 
The residue conductivity (Jres ranges 0.1-0.35 m 

g 1121241 J 

samples exhibiting electron-hole asymmetry, our analy- 
ses focus on the carrier type with the higher /xpp to avoid 
complications associated with contact doping."^ 

Although Eqs. |2] provide a good description of exist- 
ing conductivity measurements, the origin of scattering 
sources in graphene is still under debate. In addition to 
CI, ripples, resonant scatter ers and midgap states are also 
potential candidates.'^'^'^ Unlike potassium adatoms,'^ 
certain adsorbates seem to dope graphene but not de- 
grade its mobility.'^ The role of the d ielectric environ- 
ment also appears controversiaL^i^Ii^^ While Jang et al. 
find agreement with the model using ice as a top dielec- 
tric layer Ponomarcnko et al. report screening effects 
much smaller than expected from the CI model using 



liquid dielectric layers.^- Within the charged- impurity 
model, the origin of such impurities remain unclear: Ad- 
sorbates on top of graphene, charges adsorbed/ trapped 
at the graphene/Si02 interface or residing inside the sub- 
strate are all possible candidates. 
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FIG. 2: (Color online) (a) pxx(l/S) of sample A at 1.6 K and 
n = 3.89x lO^'^/cm^. The onset of the oscillations corresponds 
to filling factor v = 4:6 (arrow). The vertical bar marks 25pxx 
aXv = 22. The dashed lines correspond to Eqs.[3]with Tq — 34 
fs. (b) The corresponding Dingle plot: (5pxx/7th vs. 1/B. r, is 
extracted from the linear fit. (c) Tq(n) of samples A (squares), 
B (circles), C (up-triangles) , and D (down-triangles) on a log- 
log plot. The dashed line indicates ^/n dependence. 

We have measured the quantum scattering time Tq 
in graphene to further address the above issues. We 
determine Tq from the magnetic field dependence of 
SdH oscillations following procedures well-established in 
conventional 2DEGs.^'^ Figure |2|a) shows the magneto- 
resistance Pxx{B) of sample A at n = 3.89 x lO^'^^/cm'^ 
and T = 1.6 K. The oscillatory amplitude Sp^x can be 
described by: 



Spx 



Po 



-);ith 



2TT'^kBT/huJc 

smh{2Tr^kBT/nuJc) ' 



(3) 



Here po is the non-oscillatory background resistance, "fth 
is the thermal factor, and ujc the cyclotron frequency in 
graphene. Here to* — Ep/vp = hy/Trn/vp is the effec- 
tive mass and vp — 1 x 10^ m/s is the Fermi velocity 
in graphene. Figure |2jb) plots SpxxHth vs. l/S in a 
semi- log plot (the Dingle plot), where we extract Tq = 
34 fs from the slope of the linear fit. The correspond- 
ing 5pxx calculated from Eqs. |3]is plotted in Fig. |2|a) as 
dashed lines and exhibits excellent agreement with data. 
In each sample, the same procedure is repeated at dif- 
ferent densities for n > 1.2 x 10^^/cm^, where several 
well-developed SdH oscillations are observed before the 
onset of quantum Hall states. In some traces, a slowly 
varying background is subtracted before the determina- 
tion of 5pxx, as described in Ref.l^^l. 

Figure ^c) plots Tq{n) in samples A-D (also listed 
in Table ^ for n = 3 x 10^^/cm^). Tq{n) increases 
with increasing n in all samples and spans 25-74 fs for 
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1.2xl0iVcm^ < n < 6 X lO^Vcm^ 



These values 



correspond to a quantum level broadening T — 4.5-13 
meV and are in line with F — 20-30 meV extracted 
from the adsorption linewidth of cyclotron resonances at 
n < 1 X lO^Vcm^P 

The n-dependence of Tq in Fig. ^c) agrees qualita- 
tively with that of cr(Vg) in Fig. [T] This can be seen by 
separating the long and short-ranged components in r^.g 
using the following equations;^ 
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4| and plot the results in Figs.ji] 
Samples A and B show n- 



Due to a small p short, both tj and Tq in sample A follow 
closely the ^/n dependence expected for charged impuri- 
ties residing in the graphene plane {z = O).'^^ Large p short 
and higher p-p-^ in samples B-D (Table |l]) cause both scat- 
tering times to deviate from the -y/n dependence. How- 
ever, Tq does not correlate with in a simple relation. 
Samples B and C exhibit similar /ufe, but their Tq dif- 
fer by a factor of 2. Samples A and B show comparable 
TgS despite the significant difference in /ipE- Sample C, 
with a moderate /ipE E^nd the highest p short i exhibits the 
highest Tq. 

To understand the above observations, we calcu- 
late the ratio Tt/Tq{n) and the long-ranged component 
rj°"»/r^°"9 (n) using Eqs. ~ ~ 
(a) and (b) (also Table 
independent Tt/Tq of 2.7 and 5.1 respectively. For the 
other two samples, Tt/Tq decreases slightly with increas- 
ing n, varying from 1.9 to 1.5 in sample C and 4.3 to 
3.3 in sample D. Clearly, the angular distribution Q(0) 
in Eqs. [l] differs significantly in these samples. 

By evaluating rj°"^/Tg°"f (n), we find that such varia- 
tion may be naturally explained by varying the impurity- 
graphene distance z within the charged-impurity model. 
Theoretical calculations of short and long-ranged ra- 
tios are given in solid and dashed lines, respectively, in 
Figs. [3] (a) and (b).^ According to this model, the 
T(°"^/r^°"3(n) of sample A provides further evidence for 
the domination of CI located in the graphene plan e {z 
— 0), where a constant 2.5 is expected (dashed line) 
The Tt/Tq{n) in sample C falls between the dashed and 
solid lines, which is the result of a large short-ranged 



component. Its 



long I long 



/T^°"'3{n) exhibits an approximate 
also pointing to CI located in the 

long 1 long ■ 



constant ratio of 2.1 

graphene plane (Fig.|3[b)). rj°"^/T'°"^ in samples B and 



d^ 



D range from 5-7 and~are best described by CI located 2 
nm away from the graphene sheet (Fig. [3]^b)). In all six 



samples, we find 2 < r, 



long 



/-long 
I 'q 



< 7, corresponding to 



CI located within 2 nm of the graphene sheet . 

The knowledge of the impurity-graphene distance z is 
essential in correctly determining the CI density riimp in 
the vicinity of graphene. In the CI model, Uij^p is related 
to /ipE through /^FE= C/fiimp, where C = 5xlO^^/Vs is 
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FIG. 3: (Color online) (a) rt/r^ vs. n in samples A (squares), 
B (circles), C (up-triangles) , and D (down-triangles) along 
with calculations for CI at 2 = (dashed line) and 5-scattering 
centers (solid line) (^^). (b) Their long-ranged component 
^i^ong i^iong ^ Corresponding open symbols. Dashed 
lines are calculations for z — Q (bottom) and 2 nm (top) (f^. 

n-independcnt for z — O.'^ At a finite z, C increases with 
increasing n due to screening. Using the equations in 
j^gfgpiTj^ numerically calculate C{n, z) and estimate 
nimp by fitting a{n) to Eq. [2] The results are listed 
in Table |l] while an exemplary fitting is given in Ref.l^. 
Clearly, /ipp (or Tt) is affected by both z and riimp- For 
example, the difference in /ipE between samples A and 
D mainly stems from z, instead of riimp- In contrast, 
Tq serves as a better measure of riimp due to its weaker 
dependence on z. 

Next we briefly assess the effect of density inhomo- 
geneity 6n on the measurement of Tq . Caused by CI, dn 
measures approximately a few lO^-'^/cm^ near the Dirac 
point, and is expected to decrease with increasing n 
due to electron screening.f^S^n introduces phase smearing 
in p^^(B), effectively reducing the SdH oscillation ampli- 
tude and suppressing the value of Tq determined through 
the Dingle plot, as demonstrated in GaN 2DEGs.l^ We 
have used the intercept of the Dingle plot at l/i? = as 
a criterioiP^ to obtain Sn and the corresponding correc- 
tions to Tq. Details are given in Ref.^^Sl. We estimate 5n 
to be ~ 7 X 10^°/cm^ in sample A, and ~ 9 x 10^°/cm^ in 
sample B. These estimates are consistent with the high- 
est filling factors observed in these samples (i^ = 46 for 
sample A and v — 74 for sample B). Overall, Sn/n de- 
creases rapidly with n, in agreement with theory; but 
the magnitude is only a few percent in the density range 
we studied, which is signi ficantly smaller than the theo- 
retical predictions .fSSlSZIlll The above correction leads to 
~ 20% increase of Tq in sample A and 50% in sample B. 
The corrected TqS are given in Table [l] in parenthesis. In 
samples C and D, the corrections are smaller than the 
error bars of Tq and therefore omitted. tI°^^ /t^^^^ in 
sample B now corresponds to CI located at z = 1 nm 
instead of previously determined z = 2 nm, but the main 
picture does not changcf^^l 

Our study of Tq and tj/t^ provides critical infor- 
mation in differentiating various scattering scenarios in 
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graphene.SSEEEIl A detailed comparison to theory is only 
made for the CI model at this point, but can be ex- 
tended to other proposals as quantitative predictions be- 
come available. The diverse behavior our samples exhibit 
can all be understood very well within the CI model using 
three parameters: Uimp, z, and Pshort- We speculate that 
uncontrolled spatial variation of Si02 surface properties, 
as well as sample preparation conditions (e. g. humidity) 
may have been the primary reasons behind the observed 
differences among samples, although variations in prepa- 
ration procedures cannot be ruled out!^ 

Our results indicate that the dominant CI reside within 
2 nm of, and sometimes in the immediate vicinity of the 
graphene sheet. Primary candidates of this nature are 
charges carried by adsorbates on top of graphene and/or 
at the graphcnc/Si02 interface. The role of adsorbates is 
especially highlighted in the c urrent annealing treatment 
of suspended graphene.'^^^ In the litera ture, va rious 
approaches, including resist-free processing j-'^^ I ^^ I ^Q l UHV 
baking, ^■^ and current annealing,ISI have been used to 
remove contaminants on top of graphene without signifi- 
cant improvement to mobility. These observations collec- 
tively suggest that adsorbates on top of graphene cannot 
be the major culprit in limiting mobility at the current 
level. Instead, we speculate that charges (e. g. Na+) 
and molecular groups (e. g. OH) adsorbed/trapped on 
the Si02 surface prior to the exfoliation of graphene are 
the major source of scattering. The z = found in sam- 
ples A and C lends strong support to this hypothesis. 
Moreover, the small z observed in other samples can be 
accounted for by the existence of a spacer layer between 
graphene and Si02 (e.g. H2O). The evidence of such a 
layer is widely observed in AFM height measurements of 
graphene. In this scenario, the concentration of adsorbed 



charges (riimp), together with the thickness of the spacer 
layer (z) can account for the wide span of mobility seen in 
our samples. It may also explain why graphene on a va- 
riety of substrates displays a similar range of mobilitjESl 
since the bulk properties of these substrates are less rel- 
evant here. 

Finally we note that the above determined z can be 
expanded to represent an average impurity-graphene dis- 
tance. Using this concept, we consider the contribution 
of uniformly distributed charges within the bulk of the 
Si02 substrate. Our simulations show that oxide charges 
in commercially available Si02 are unlikely to be a major 
source of scattering at the current mobility level;-- 

In conclusion, we have systematically studied the 
quantum and transport scattering times in graphene. 
Our data will prove useful in critical examinations of ex- 
isting scattering scenarios. Within the CI model, the 
ratio of Tt/r^ indicates that charged impurities residing 
within 2 nm of the graphene sheet are the main sources of 
scattering in graphene. Such information provides impor- 
tant guidance to the effort of improving carrier mobility 
in graphene. 
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